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Following a previous study of the one-loop factorization of the nonsinglet g 2 structure function 
of the nucleon, we present in this paper the next-to-leading order coefficient functions in the 
singlet sector. To obtain the result, the partonic processes of virtual Compton scattering off two 
and three "on-shell" gluons are calculated. A key step in achieving the correct factorization is 
to separate the correct twist-two contribution. The Burkardt-Cottingham sum rule is nominally 
satisfied at this order. 



In lepton-nucleon deep- inelastic scattering with the electromagnetic current, two spin-depen- 
dent structure functions of the nucleon can be studied: 51,2(^5, Q 2 )- gi(xB,Q 2 ) is closely related 
to the spin structure of the nucleon and has been investigated extensively in the last decade. 
g2{xB,Q 2 ) is present in processes involving transversely polarized nucleons and is a three-twist 
structure function in the sense that it contributes to physical observables at order 1/Q M. The 
experimental measurements of the g 2 structure function have been done by several collaborations 
@. There is a long history in debating the physics involved in g 2 [xBiQ 2 \ As more results 
from quantum chromodynamics (QCD) are obtained and understood, it is generally accepted that 
g 2 {x b,Q 2 ) probes quark and gluon correlations in the nucleon which cannot be accessed through 
Feynman-type incoherent parton scattering. 

At leading-order, g 2 (xB, Q 2 ) can be expressed in terms of a simple parton distribution Aqx(x), 

g T (x B , Q 2 ) = (gi + g 2 )(xB, Q 2 ) = \Y. e l (MaAxB, Q 2 ) + Ag aT (x B , Q 2 )) , (1) 

where 

Agr(x) = ±J ^(PS ± \mi±l^n)\PS ± ) . (2) 

This result has lead to many incorrect interpretations of the gi physics. When the leading- 
logarithmic corrections are studied, it is found that Aqr(x, Q 2 ) mixes with other more complicated 
distributions under scale evolution ||. In fact, Ag^(x, Q 2 ) is a special moment of a general class of 
parton distributions involving two light-cone variables [|], ||. When the scale changes, only those 
general distributions evolve autonomously. 

Thus a general factorization formula for g 2 is much more involved than the leading-order 
result shows. It should contain the generalized two-variable distributions, Ki(x,y). Indeed, we 
shall write in general 

9t(xb,Q 2 ) =J2 f ——{Ci(—,—\a a ]Ki(x,y) + (x B ^-XB)\, (3) 
i J-i x y { \ x y J J 

where Cj are the coefficient functions. In a previous paper , we studied the factorization of the 
nonsinglet part at the one-loop order where there were two distributions, K\ a ^ associated with 
each nonsinglet component a, and their one-loop coefficient functions were obtained for the first 
time. 

In this paper, we study the one-loop factorization of the singlet part of gr{xB, Q 2 )- The subject 
has previously been considered in Ref. [0, a comparison of the results will be made in the end of 
the paper. Contrary to the previous conclusion, the result of this paper represents a local operator 
product expansion. Throughout, we use the kinematics defined in ||. We first define the singlet 
quark distributions, 

Kix(x,y) = Y / K ia (x,y) (4) 
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Figure 1: Feynman diagrams contributing to two quark and one gluon intermediate states. The 
intermediate longitudinal gluon is onshell, so a special propagator derived from the equation of 
motion has been employed. There are a total of six diagrams. In dimensional regularization, the 
sum is zero. 

where the sum is over all quark flavors. Their one-loop coefficient functions are 

C?U^M=?y5(x-x B ), (5) 
\ x y J 

where e 2 = J2i e l/^f is the mean square quark charge and Nf is the number of active quark 
flavors. Using the relation 

Ag T (x) = - fdy{K l {x,y)+K 2 (x,y)} , (6) 
x J 

it is easy to see that the result is consistent with Eq. (1). 

The general strategy of higher-twist factorization beyond the leading orders has been presented 
in || and will not be repeated here. For the singlet factorization, we need to consider four classes 
of parton intermediate states: two quarks, two quarks-one gluon, two gluons, and three gluons. In 
the cases of two-quark and two-quark-one-gluon states, many of the relevant diagrams have been 
considered in 0; we will not repeat that analysis here. In addition we must take into account the 
explicit singlet diagrams shown in Fig. 1. A detailed calculation shows that their contributions 
cancel. Thus the one-loop coefficient functions of (x, y) remain the same as those of the 
non-singlet sector. In what follows, we focus entirely on two- and three-gluon intermediate states. 

The diagrams involving two-gluon intermediate states are shown in Fig. 2. To isolate the 
0{1/Q) contribution, the incoming gluon is given a transverse momentum which is expanded to 
first order in the internal propagators. Our calculation shows that only the correlation function 

n gB (x,y) = 8{x -y)j (PS\A^O)id a A a a (Xn)\PS) (7) 

contributes, where a and /i take only transverse values. 

The perturbative diagrams corresponding to three gluon intermediate states are shown in Fig. 
3. The relevant correlation function is 

I% B (*,V) = J ^^e lA V^(P5|(-z)^/ afec ^(0)A a (^)^(An)|PS) . (8) 
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Figure 2: Two-gluon contributions to gr- The cross <g> represents one transverse momentum 
operator insertion. There are a total of eight diagrams. 





Figure 3: Feynman diagrams for three-gluon Compton scattering. There are a total of 24 
diagrams represented in (a). The other 8, in figure (b), have the form of those in Fig. 1, and 
vanish for the same reason. 

All fields and couplings in the above expressions are bare. Gauge invariance demands that the 
final result depends on the combination 

= J ^^e iA V^)(P5|^(0)^ Qa5 (/in)^(An)|P5) (9) 

of distributions, where D^ ab = <5 a 6<9 M — gf ac bA c . In addition, gauge invariance demands that 
physical observables depend on moments of 

xyr» B (x,y) = /^^e a V^)(P5|F^(0)^ Qa6 (HF^(An)|P5) 

= K gB (x,y)te a ^ Pa n p S, . (10) 

With appropriate insertions of light-cone gauge links, which can be generated by summing over 
states with additional longitudinally-polarized gluons, K gB is gauge invariant. Our goal is to 
calculate its one-loop coefficient function. 

The Compton amplitude for virtual photon scattering on a transversely polarized nucleon can 
be written as (e 0123 = +1), 

T v> = -ie^q a S ±f} -S T (x B , Q 2 ) . (11) 
v 
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In QCD, we can write the singlet part as 

" l dx dy 



S^(x B: Q< 
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-xb 



where the M's are perturbation series in a s and have infrared poles. Once again, the tree and 
one-loop level expressions for Mj£ are the same as the nonsinglet case in Ref. ||. 

The amplitude M g starts at the one-loop level. To simplify the expression, we assume \x B \ > 1 
so that the amplitude is purely real. 

Since all Feynman diagrams are computed in bare perturbation theory, the result depends 
on the bare coupling gs- We replace it with the renormalized coupling in the MS scheme; the 
difference appears only at higher orders in a s (fi 2 ). A straightforward calculation yields 

a s (^)T F ftfY' 2 



MP(x,y) 



An 



,Q 2 



1 

2xy 



x 



((y - 2x)(l -x) + 3y(y - x) - 6 | (y - x)j log (1 - x) 



x 



+ [(Ax - 3y)(l - y) + Ax(y - x) - 8- (y - x) \ log(l - y) 



+ 



(x + y)(2x-3y) 



x-y 



-3[y + 2 



(y - xf 



(1 - (x - y)) - xy log (1 - (x - y)) 



x) logfl — x) 



X 



+ \^3(Ax - 3y)(l -y) + Ux(y -x)-20-(y- x)\ log(l - y) 

6x 2 + xy - 9y 2 . . . 

H (1 - (ar - y>) log (1 - [x - y)) 

x-y 

~2 V V ~ 2X ^ 1 ~ x "> + 3y ( y ~ x ^ ~ 6 | ( y ~ x log2 ^ - x ^ 
-~ ((Ax - 3y)(l -y) + Ax(y -x)-8^(y-x)J log 2 (1 - y) 
_l_f(x + y)(2x-3y) 



x-y 



(1- (x- y)) - xy log (l-(x- y)) 



(13) 



where e = 4 — d, = Ane'"" 3 ^ 2 , and 7^ is the Euler constant. We have also introduced the 
generator normalization Tp = 1/2. 

Now we want to show that S^(xb,Q 2 ) is factorizable at the one-loop level, i.e., the infrared 
poles 1/e in M g match the ultraviolet poles in K 9 b- To this end, we use the infrared poles in M g 



A 



to generate a scale evolution equation for the renormalized parton distributions, 



d log /i 2 



1 / 2 
dxdy — 

-1 \X[Xb — X) 



a 



j (K^(x, y, /i 2 ) + K 2S (x, y, /i 2 )) 
dxdy 



8n 7-i x 2 y 2 (x — y) 
({y - 2x)(x B - x) + 3y(y - x) - 6 (y - x)^j log N 



x 



x B 



+ [(Ax - 3y) (x B -y) + Ax(y - x) - 8 (y - x) log 1 

V y \ x B 



+ 



(x + y)(2x - 3y) \ / ./■ - // 
(x B - (x - y)) - rry log ( 1 



x- y 



K g (x,y) + 



(14) 



where the ellipses denote the homogeneous part of the evolution. This contribution is given 
explicitly in Ref. |J. In order to compare this result to the known twist-three evolution ||], we 
must first remove the twist-two contribution. Consider the twist-two operator 



(15) 



whose indices have been symmetrized and whose traces have been removed. Its matrix elements 
in the nucleon state are given by 



(PS\9^ + T n+1 \PS) = 2a n P { ^ ■ ..p^s" 



n + 1 



(16) 



It is the _!_+•••+ component of 9 that appears in this transverse process : 



g(i+-) 

'n+l 



c % 3 



n + l 



(n + l)iA\id + ) n id k A k + iA J (id + ) n - L D a F 



n-l 



+ f abc A j a {id + ) n ^gA k b id + A k c + f abc Ai(id + ) m gA k (id + ) n - m A k 



m=0 



(17) 



Here, latin indices are understood to take values in the transverse dimensions. Expanding Eq. QT3| ) 
in the large xb limit, one arrives at evolution equations for the moments of the parton distributions. 
Removing the twist-two part of these equations via Eq. ([17]) , one can obtain the evolution of the 
twist-three operators. A detailed check shows that our result is identical to that found in Ref. j|] 
obtained by studying the ultraviolet divergences present in the twist-three operators. We note here 
that this separation provides a new homogeneous term in the evolution of the singlet twist-three 
quark operators. In the absence of this term, the diagonal evolution of these operators would be 
identical to that in the nonsinglet sector since the contributions displayed in Fig. 1 vanish. 

The final step of the calculation is to take the imaginary part of the factorized S^(xb, Q 2 ) to 
get a factorized expression for the structure function g^(xB,Q 2 ) in the physical region x B < 1- 
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We find the coefficient function 



£,(1) / x B x b 
9 V x ' V 



a s (fi 2 )T F 1 
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x-y 
x B 



1 



where 0(x) is the step-function, which appears in Eq. along with a factor of l/(y — x). 

To check the Burkhardt-Cottingham sum rule ||, we integrate g^{xs, Q 2 ) over xb- Assuming 
the integration over x and y can be interchanged with that of xb, one obtains 



dx B gT{xB,Q 2 



h 2 1 7 



_ a s (Q 2 ) 
2 F 2tt 



(PS'lE^mTsV'il^s) 



(19) 



where 75 has been defined in the 't Hooft-Veltman scheme. The coefficient 7/2 reduces to 3/2 if 
we define 75 so that the nonsinglet axial current is conserved. Compared with the factorization 
formula for gi(xB,Q 2 ) H, we have the Burkhardt-Cottingham sum rule at one loop; 



/ dx B gf(x B ,Q 2 ) = . 
Jo 



(20) 



If the order of integration cannot be interchanged because of the singular behavior of the parton 
distributions at small x and y, the above sum rule may be violated. Indeed, some small xb studies 
indicate such singular behavior [ |T0|j . 

Finally, we consider the next-to-leading order correction to the singlet part of the x 2 moment 
of gx{xB, Q 2 )- In the leading order, it is well known: 



jf 1 dx x 2 ^(x, Q 2 ) = h 2 Qa 2S (Q 2 ) + d 2E (Q 2 



f21^ 



where a 2 is the second moment of the <?i(x, Q 2 ) structure function and g? 2 is a twist-three matrix 



element [11|. Using the coefficient functions found above, we have 



dxx g T (x, 



-e 2 



'a 2S (g 2 



1 + 



4tt 12 



a 29 {Q 2 ) a s (Q 2 ) 5 



Aix 



N f T, 



f i F 



+d 2 ^Q 2 



1 + 



a s {Q 2 ) /27 



47T 



T c A - 2 i CF+ ^N fT ; 



.(22) 
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Using the next-to-leading result for gis(x,Q 2 ) [[|, we find 
J'dxx 2 ($(x, Q 2 ) - ^gf(x, Q 2 )) 

= h?^ m |l + (?L Ca - fc F + f N f T F ) } . (23) 

Notice that the combination of Qt and g\ relevant to c?2s receives no radiative correction, as in 
the nonsinglet sector. A detailed calculation shows that analogous results are valid for all higher 
moments as well []12] . This implies that the tree relation between gi(x,Q 2 ) and the twist- 2 part 



of g 2 (x, Q 2 ) is respected at one- loop order. These results and their implications will be presented 
in a future communication [I2"|. 



The significance of the present result is as follows. In the leading order analysis of g 2 , one just 
needs the leading- logarithmic evolution of Ki(x,y) which is now well known j|, |5j, |13|| , including 



its large N c behavior . In the next-to- leading order, one needs to know the coefficient functions 
and the two-loop evolution of Ki(x,y). The former is now complete with Ref. and the present 
paper. The latter has not yet been calculated, but in generally its effort is not as important as 
the coefficient function we calculate here. 

Finally, let us add a few remarks on the comparison of our findings with a previous calculation 



of Ref. J7|. As has been noted earlier [[Uj the result of that paper is not complete since the 
contributions of twist-three two-gluon operators, i.e. the diagrams on Fig. |2|, were not accounted 
for. Since the effect of two-gluon graphs on the final answer is reduced to a redefinition of the 
correlation function (schematically) from {A^AjjA^) to (A^D^A^ ) without affecting the three- 
gluon coefficient function, we can directly compare both results. Before we do it, we observe that 
the basis of three-gluon functions used in |fTE, [^| is redundant, namely, in the decomposition 



(A^AyAp) oc N(x 1 ,x 2 )g p L p e u - +a s^ + N(xi, x 2 )s^e ILV ^ + + (24) 

obviously, the Lorentz structure in front of N can be expressed in terms of the first line using an 
identity 

Therefore, in final equations of we have to put N = to be consistent with the present analysis. 
Moreover, the result in Ref. [|16| must be also modified correspondingly. 



Next, to have a correspondence with the correlation functions used here and in Ref. H, we 
have identify the momentum fractions of the gluon lines as follows: x = X\ and y = X\ — x 2 - For 
the latter convenience we introduce as well a new coefficient function 

£^(x,y;x B )^^-C^(^M. (26) 
y y — x y \ x y J 
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Since the function N(x\,x 2 ) is symmetric w.r.t. interchange of its arguments only symmetric part 
of £^'(xi,Xi — x 2 ;xb) will survive. Thus, a symmetrization leads to the equation 

£^\xi,X! -x 2 ;x B ) + £^(x 2 ,x 2 -x x ;x B ) = ^E 2 (x 1 ,x 2 ;x B ), (27) 

07T 

where E 2 is the coefficient function from Ref. 
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